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Abstract. Abstract. A well-posed initial-boundary value problem is for- 
mulated for the model problem of the vector wave equation subject to the 
divergence-free constraint. Existence, uniqueness and stability of the solution 
is proved by reduction to a system evolving the constraint quantity stati- 
cally, i.e., the second time derivative of the constraint quantity is zero. A 
new set of radiation-controlling constraint-preserving boundary conditions is 
constructed for the free evolution problem. Comparison between the new con- 
ditions and the standard constraint-preserving boundary conditions is made 
using the Fourier-Laplace analysis and the power series decomposition in time. 
The new boundary conditions satisfy the Kreiss condition and are free from 
the ill-posed modes growing polynomially in time. 



1. Introduction 

When formulated for the purpose of numerical integration the equations of gen- 
eral relativity are customary divided into the two subsets: the subset that contains 
both time and space derivatives, called the evolution subset, and the subset that 
has spatial derivatives only, called the constraint subset. Usually, the evolution 
subset is the one that is being solved, while the constraint subset is made automat- 
ically compatible by a choice of initial and boundary data. Exact compatibility, 
however, can never be achieved in numerical calculations, which are the primary 
source of solutions to Einstein's equations. To predict accuracy of the solution, a 
new formalism is needed to describe systems of constrained evolution. 

In fact, it has been a challenge in numerical relativity to keep the constraint 
subset of equations from exponential growth [SJ [12l [21] [30] , when evolution equa- 
tions are solved exclusively. Methods were introduced to enforce compatibility of 
the solution with the constraint equations by adding differential terms to the evolu- 
tion equations so as to minimize growth of the constraints (e.g., [3T]), periodically 
resolving constraint equations, when they become too large, (e.g., |14j). and con- 
trolling the constraint perturbations from entering the domain through the timelike 

boundaries (e.g., [HHE1E1H1ES1E1E1I25])- The method proposed in this paper 
is related to the last group. We explore the idea of [10] and study the closely related 
system that propagates the constraint quantities statically. The results are used to 
formulate a well-posed initial-boundary value problem for the original constrained 
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evolution system. On the next step, the constrained problem is used to define a 
well-posed initial-boundary value problem for the free evolution. 
We consider a model problem of the vector wave equation 

(1) dfu t = WdjUi, 
subject to the differential constraint 

(2) d l Ui = 0. 

The choice of the model problem is motivated by re-formulations of Einstein's field 
equations (e.g., [3l [34J [4l [3ll [24l [TTJ [12l [26] ) in which second order in space equations 
are coupled to first order differential constraints (if necessary, ([T]) can be reduced 
to first order in time by introducing tt{ — dtUf, all results can be repeated for 
the pair (u% , ~Ki ) without much complication) . First order symmetric hyperbolic 
reductions of UJ has been used as a source of insight in numerical relativity in 
the past [6j [9] [211 HZ] and methods for constructing constraint-preserving boundary 
conditions were investigated using the same model problem. In this work we prove 
that (TTJ), ([2]), taken as a constrained evolution problem, is related very closely to the 
equation 

d^Ui = d^djUi — did^Uj 

which evolves the constraint ^ statically. Based on this observation, we formulate 
an initial-boundary value problem for the constrained evolution system (UJ), ([2]), 
and prove existence, uniqueness, and stability of its solution. 

The first application of the new approach is the new set of radiative constraint- 
preserving boundary conditions (see e.g., [21] [16j HZl E 9 , 18, 28 ) for the free evolu- 
tion problem (TTJ). In the simplest case, when homogeneous conditions are specified 
on the incoming waves, the boundary conditions take the form 

(3) d t u A + d n UA - dAU n = 0, d t u n + d n u n = 0, 

where d n stands for the normal, and 8a for the tangential derivatives at the bound- 
ary, and ti„, ua are projections of Ui on the outward pointing unit normal and the 
tangential directions, respectively. 

In |27j , an example of constraint-preserving radiation-controlling boundary con- 
ditions is constructed for a symmetric hyperbolic reduction of the fat Maxwell 
system. The boundary conditions ([3]) overlap with the boundary conditions of |27j 
in two, but differ significantly in one equation. Following the standard approach of 
the subsidiary evolution system for constraint quantities, the radiative boundary 
condition enforcing the constraint in [27] takes the form of a first order differen- 
tial condition (containing both spatial and temporal derivatives) for the first order 
system. Here, however, the new approach allows us to write all three conditions 
as first order conditions (but still including tangential derivatives) for the second 
order system, and existence is proved without any additional gauge assumptions. 

The problem (jTJ) , ([2]) has been considered in [18] in the case of a four-dimensional 
vector it, e.g., u = (ut,u x ,u y ,u z ). As in [27], the intuition derived from Maxwell's 
equations has been used to formulate a set of Sommerfeld-type boundary conditions 
that are constraint-preserving. In fact, three of the four conditions of |18j coincide 
with ([3]) if the component Ut is set to zero. The fourth condition, however, is closely 
related to the homogeneous Dirichlet data for the constraint quantity. It does not 
vanish for u t = 0, and, therefore, the boundary conditions of |18) do not reduce to 
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Several authors (e.g., |16[ I26j) had experimented with boundary conditions mo- 
tivated by physics rather than the equation analysis. Surprisingly, the stability of 
the calculations improved when physical data was employed. Their results moti- 
vated us to search for the mathematical framework that could explain such effect, 
and possibly describe the accurate boundary conditions for constrained evolution 
systems. 

The idea of reducing Einstein's equations to a form evolving the constraint quan- 
tities statically was first used in [TD] - By exploiting gauge freedom, the authors 
re-write the tetrad formulation of Einstein's equations as to have the constraint 
quantities propagating along the boundary (but not inside the domain) , thus mak- 
ing any boundary conditions that are well-posed for the problem automatically 
compatible with preservation of the constraints. This approach allowed the au- 
thors to formulate a set of maximally dissipative boundary conditions for which 
well-posedness can be proved in the full nonlinear case using standard theorems for 
symmetric hyperbolic systems. 

In this work we use constraint equations — rather than gauge freedom — to for- 
mulate the static constraint evolution reduction, as described in Sections 2 and 3. 
In Section 4 we use standard theorems for first order symmetric hyperbolic systems 
to prove the existence of the solution to the static constraint evolution problem. 
In Section 5 we formulate a well-posed initial boundary value problem for the con- 
strained evolution system (UJ), ([2]). The existence, uniqueness, and stability of the 
solution follows from the theorems for the static constraint evolution problem. We 
also formulate a set of constraint-preserving radiative-controlling boundary condi- 
tions for the free evolution problem (fT]) and prove that this set is well-posed. Notice 
that the norm in which the well-posedness is proved is built on the antisymmet- 
ric part of the gradient and the divergence of the solution. However we believe 
that this is the natural norm for the problem, and any "extra" estimates on the 
gradient require additional assumptions, like more restrictive boundary conditions 
exemplified in Section 6, or a special gauge choice similar to one in [27] . 

Section 7 is devoted to the Laplace-Fourier analysis of the free evolution problem. 
We prove that the new boundary conditions ([3]) satisfy the Kreiss condition, thus 
the free evolution problem is stable for zero initial data. (Notice that the stability 
of the system can be studied for arbitrary initial data using techniques of fL7\ ITS] 
if the Kreiss condition is satisfied) . We find that the Neumann-Dirichlet constrain- 
preserving conditions obtained by the standard method do not meet the Kreiss 
condition. Also, the radiative differential conditions following from the standard 
approach were found not to meet the Kreiss condition. Finally in Section 8 we 
study stability of the system in the case of non-vanishing initial data using the 
framework developed in |28] to quantify instabilities that grow polynomially in 
time. It was found that both the new conditions ((3]) and the radiative differential 
condition similar to one in [27] do not expose such instabilities, while the standard 
Neumann-Dirichlet data does. 



Notations 

We denote by a bounded convex subset of R 3 with boundary consisting of parts 
of planes. The polyhedron domains are relatively easy to implement numerically, 
and therefore represent an important example from the practical point of view. 
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Vector m in most cases denotes the outward pointing unit normal vector to d£l. 
Vectors rm and U are tangential to dQ and complement to an orthonormal triple. 

We use notations of general relativity, in which the repeated indices denote 
summation, small roman indices run from 1 to 3, the indices are raised and lowered 
with the flat metric £y = {1 if i = j, if i ^ j}. Later we use parenthesised and 
bracketed indices to denote the symmetric and anti-symmetric parts of an array 
correspondingly, e.g., = (Wij+Wji)/2, W[ tj ] = (i% -Wji)/2, Wij = w^+w^j]. 

The symbol di denotes the gradient, however, d n is reserved for the directional 
derivative in m. We will use capital roman indices to indicate the tangential com- 
ponents of the fields, for example, &a stands for derivatives in tangential directions 
to dfl (or, if rii is not associated with <9f2, Da will denote derivatives in directions 
normal to ni). Similar notations are used for vectors: v n — n 3 Vj, va — (Sj — riin J )vj. 

2. Subordinate constraint evolution and constraint preservation 

Our goal is to formulate a well-posed initial-boundary value problem for the 
constrained evolution system ([I]), |2]). We introduce the constraint quantity 

C = d l Ui . 

Lemma 1. If u i is a solution of then the constraint quantity C obeys the wave 
equation, 

(4) d 2 t C = &diC. 

Proof. The lemma follows by taking the divergence of Eq. ([1]) and commuting 
derivatives. □ 

Corollary 1. Consider the problem of solving Eq. fl]] in R 3 . Let Ui be a solution 
to (Qp, corresponding to constraint- compatible initial data, i.e., C(0) = d i Ui(0) = 
and dtC(0) — d l dtUi{§) = 0, and Ui is regular to guarantee existence of the spatial 
integrals of C and its derivatives. Then, C = 0. 

Proof. According to Lemma [lj the constraint quantity C is a solution to (j4|). Mul- 
tiply (jlj by dtC and integrate over R 3 . After integration by parts, 

/ [dfC){d t C)dx = - j (diCWdt^dx. 

By commuting derivatives and re-grouping terms, we obtain 

\dt[\\d t C\\l 2m + ||^<7||! 2(a3) ] = 0. 

Therefore, ||5tC||^ 2 ^ R3 -j + ||^C||^ 2 ^ R3 ^ is constant in time. Since it is zero initially, 
it remains zero for all times, which in turn gives that C = 0. □ 

Remark. Note that Corollary Q] does not hold if M 3 is replaced with a bounded 
region. Indeed, if C ^ at some point of the boundary, in view of ((3]), the constraint 
perturbation will flow inside the domain regardless of the compatibility of the initial 
data. 

In general to ensure C = the boundary data must be given consistently with 
the evolution of the constraint quantity. The data that guarantees that C remains 
zero for all times is called constraint-preserving. 
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3. Static constraint evolution reduction 
In domain fi, consider the equation 

(5) $1* = WdjUi ~ diC 
subject to constraint ([2]). 

Lemma 2. Evolution of the constraint quantity C in Eq. {5[) is static. That is, let 
Ui be any (sufficiently smooth) solution to then 

(6) d 2 t C = 0. 

Proof. The proof follows by substituting the definition of C, commuting derivatives, 
and using ([5]). □ 

Remark. One may notice that in ([5|), we recover the second order reduction of 
Maxwell's equations in the Coulomb gauge. Indeed, by substituting the definition 
of C, we rewrite (0 as 

d 2 Ui = 2d 1 duUa = (curl curlw)i, 

where (curlu); is the usual curl operator defined by (curlu)i = e/ k djUk', here e^-fc 
is the unit totally antisymmetric tensor, and we have made use of the identity 
£ P jk£ p - 2d [j d k]- 

It is known that the divergence constraint is preserved trivially for Maxwell's 
equations. 

Using notations of [33], we introduce the linear algebraic operator L^f q = 2<5|^ 
acting on matrix fields in f2 C R 3 , and rewrite Eq. ([5|) as 

(7) flfo = d'L^dpUg. 

Lemma 3. The operator Ejf q is symmetric in the sense of the scalar product 

(4>ij,ipij) ■= / ^ ij ip ij dx= / fajijjabS^S^dx. 
Jn Jn 

Let vectors ni, mi, and U constitute an orthonormal triple in R 3 . The eigenvalues 
and eigenvectors of L-^ q are given by 

2: y/2 V2 

: V2 f^jTOj), V2n(jli), v^m^jl^, HjUi, rrijmi, Ijli. 

Proof. The symmetry property follows by direct substitution and rearrangement of 
dummy indices: 

(^,L/V P9 >= / ^Lj q iP Pq S a iS M = f fait* 6* -5^5^ pq 
Jn Jn 

Jn 

Note also that in the standard basis for 3x3 matrices, el = riifij, e2 = niJUj, 
e3 = nilj, . . . , e9 = 14 j, the operator L-f q takes the form of a 9 x 9 symmetric 
matrix. □ 
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Remark. Spectral structure of Ljf q suggests that Eq. ((7j) does not allow control of 
the components of d p u q that are spanned by the zero eigenvectors, since Ljf anni- 
hilates them. However, this conclusion is often only partly true. An energy estimate 
on the entire gradient can be established using methods described in Section 6. 

4. First order reduction and boundary conditions 

Well-posedness of ([7j) , ([2]) is proved by the first order symmetric hyperbolic re- 
duction. We introduce the components of djUi spanned by the nonzero eigenvectors 
of Lj? q as the new variables, namely, ipji = 2d[jU^. Then Eq. (|7|) is reduced to 





d t Ui 


= 7Ti, 




dm 


= Plpji, 


(8) 







The initial data for the new variables is given by 



(9) 7r < (0)=a t « i (Q), ik<(0) = 20^(0). 

Remark. Equation ipji = 2dyUji, the definition of variable ipji, formally has to be 
added to © as an extra (artificial) constraint equation to guarantee the equivalence 
of © and ([7]). However, this constraint follows from the last equation of © by 
integration in time. Because of the simple behaviour, we postpone it in the analysis. 

Remark. A different first order decomposition, corresponding to ipji = djUi, iri — 
d t Ui, can be proposed for (J7J). The resulting system is called the fat Maxwell system 
and is only weakly hyperbolic. A detailed study of this system can be found in [27j . 
However, Lemma 2 predicts that the symmetric part of the variable ipji is not 
"seen" by the right side of ((?]). This results in the incomplete set of characteristic 
variables for the fat Maxwell system. It is possible, however, to obtain a strongly 
hyperbolic system by adding a suitable combination of constraints to the evolution 
equations. 

Let us introduce the scalar product for vectors fields, (v,i,Vi) — f^ifvidx, and 
the norms ||«||f, a ( n ) = IMIt^n) = {i'ij^ij)- B Y contracting the last 

two equations of ((5J) with 7r l and ip Jl correspondingly, integrating by parts, and 
re-grouping terms, we have 

(10) $[2|Mli a( n) + IMli a (n)l = 4 / nyn^da. 

The boundary conditions for system (jSJ) are obtained from the theory of sym- 
metric hyperbolic systems. Let m be the outward unit normal to dQ. Then the 
characteristic speeds and fields in the direction m are 

: 7T„, ipAB, 

±1 : = ir A ± ip nA . 

On each face of the boundary <9°,, we impose the maximal non-negative boundary 
conditions (in the sense of [OS E3 H3 [22] , see also in [2TT [27]^. 

(11) w\ = a A Wg +g A , 

where the smooth matrix a\ is defined on 90, and satisfies ||a|| = sup ^[^nr-^ < 1- 
Since one may think of the characteristic variables as the coefficients of the zero, 
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incoming, and outgoing modes correspondingly, condition relates the amount 
of incoming radiation to the outgoing. 

Let us consider the case of homogeneous data first. If qa = 0, then (fTTj) implies 
that the right side of (p~0|) is non-positive. Indeed, introducing the notation \v\ 2 = 
vav a 7 we have 

(12) 4ir A ip nA = \tta + ^ nA \ 2 - \tt A - < (\\af - l)\ W -\ 2 < 0. 

Let the differential operator L be defined by Eq. ([5]), i.e., 

L(tt u ipji) = (d t TTi - cPipji^tipji - 29y7r i ]). 

Let 

n L = e L 2 (n x [o,T]) | LfaM e L 2 (n x [o,r]), 

||(7Tj,^i)l|-Hi = IK^^jjOIUaCHxIO.T]) + ll-^^ii V'j^lU^nxIO.T])}) 

where ipji)\\L 2 {nx[o,T]) = |K;|U 2 (nx[o,T]) + ||V'ji||i2(nx[o,T])- It was shown in 
[25] (see [25], Theorem 2) that the integration by parts formula holds for (m, ipji) G 
Hl, which implies that the boundary integral in ([Tu]) is well-defined. Moreover, 
one can apply the result of [55] , to formulate the following 

Theorem 1. A unique solution G TLl exists satisfying {#[), £9[) and ill]) 

(with gA = 0j. In addition, (7^,7/;.^) G C((0,T) : L2(^)) and satisfies the estimate 

sup [2||7r(t)|| £a(n) + \\mh 2 (n)} < [2|k(0)|U 2(n ) + ||^(0)|U a(n) ]. 

0<t<T 

Remark. The original theorem in [25] uses the assumption that the spatial boundary 
is C 1 while here dft is only piece- wise smooth. However, in [25] it is suggested how 
to generalize the result to a boundary with corners. In general, both initial and 
boundary data have to be compatible at the corners to maintain smoothness of the 
solution. 

Remark. A similar statement can be proved using methods of JT5J dHl HO] [53] . 

In the case of inhomogeneous boundary data (gA(x,t) 7^ 0) we will assume as 
in [25] that the problem is reducible to the homogeneous case by subtracting an 
H 1 function from the solution. For that, we require that the function gA (x,t) G 
H 1 ^ 2 (dfl x [0, T]) and satisfies the necessary compatibility conditions at corners so 
as to ensure that a vector field §i G i7 1 (f2 x [0,T]) exists such that gA\on = 9a/^- 
Moreover, we assume that a vector field g\ G _ff 1 (il x [0, T]) exists such that (a) 
g A 9A\dn = 0; (b) at each point of boundary, \\cja\\ — ||<?a||/2 (here |j • |j is referencing 
to the length of the two- vector). From these two assumptions, it is not difficult to 
deduce that 

n J SjiAg l = (±)gA = (±)5a/2, on dQ. 

We will assume that coefficients of gi are chosen in such a way as to have nPejiAg = 
9A = gA/2. 

It is straightforward to verify that if fields gi and gi with the above properties 
exist, the problem ([8]), ([9]), and rtTTj) can be reduced to the homogeneous case (with 
a non-trivial forcing term in the evolution equations) by introducing the variables 

TTi = TTj - g it •ijjji = ipji - Sjug 1 . 
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Remark. In Sections 7 and 8, methods are presented to study properties of the 
solution in the case of inhomogeneous data g A in pip not invoking the above 
reduction. Instead, the techniques of the Laplace-Fourier transform are used for 
the half-plane reduction of the problem [13j [17j [18] [28] . 

The result in [55] implies that boundary integral in (flT)]) is defined. However, in 
(fTTjl holds, we can strengthen the result slightly by showing that, for (7^, ipji) € Hl, 
the traces w\ = tt a + ip nA and w A = tt A — ipnA are L 2 {d^l). 

Lemma 4. Let a vector field TTi, and an anti- symmetric matrix field ipji be such 
that TTi € L 2 (fl), duTt{i £ L 2 (fl), and similarly, ip^ £ L 2 (fl), d^ipji G L 2 (£l). If, in 
addition, 7Tj and ipji satisfy ill]) with g A = 0, then 

\\n A - Vwi||| 2 ( a n) < y— jj^p[IMlL(o) + H^Wl^n) + IMlL(n) + ll^^lll^n)] 

: = 1 _ || a || 2 [Ikllg^curl^) + IMIff^div^)]- 

Proof. The following identity can be verified by integration by parts: 

d\jTTi]i) ji dx + / d j ij)jiTT i dx= I Tr A ijj nA da. 
', Jn Jan 



Recalling ([12]), 

/ 7r A r A da >(l-||a|r)||^-^A|| 2 L2(a o)- 
Jdn 

The statement of the lemma follows by applying the Schwartz inequality to the left 
side of the integral identity. □ 

Adjusting notations of [25 to the variables of system ([5]), one can formulate the 
following theorem 

Theorem 2. Let conditions Ml)) be imposed on the boundary dfl x (0, T) and let 
the vector function g A (x,t) £ H 1 ^ 2 (dil x [0,T])) be such that the fields gi and gi 
exist with the above properties. Then there exists a unique solution (ni,ipji) € Hl 
to {3[) 7 @). I n addition, (iTi,ipji) € C((0,T) : L 2 (Q)) and satisfies the estimate 

sup [2||7r(t)|U s( n) + ||V(t)lk(n)] < c I ||flA||ffV2(an)*+[2||7r(0)|| L2(n) +||V(0)|| L2( n)], 

0<t<T Jo 

where the constant c is independent of TTi , ipji . 

Remark. A similar statement for ([8]), ([9]), ([TT]) can be found in [23] . 

5. Constraint-preserving boundary conditions for the second order 

equations 

When written in terms of the original variables of Eq. ([7]), the characteristic 
speeds and variables read 

0: <9 t it„, d[ A u B] , 
±1 : w A = d t u A ± 2d[ n u A] . 
The maximally dissipative boundary conditions (jllj) take the form 
(13) (d t u A + 2d [n u A] ) = a A (d t u B - 2d [n u B] ) + g A . 

Using Theorem [2] one can formulate the following 
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Theorem 3. There exists a unique solution to the problem Fty, \13\) . Moreover, 
the solution obeys the energy estimate: 

sup [Wdtu^L^n) +2||%ti i] || i3( n)] 

0<t<T 

(14) <c / ||9A|lHV 2(a n)^ + [||^(Q)IU 2 (n)+2||%^ ] (0)|U 2( n)]. 

Jo 

We use well-posedness of ([7]), (|13|) to establish the following result. 

Theorem 4. The constrained evolution problem (Q]) ; fl|], provided with boundary 
conditions is well-posed, i.e., for any initial data compatible with the constraint 
(in the sense of Corollary\^j), a unique solution to (Tip, I113\) exists and satisfies 
the estimate of Theorem^ 

Proof. To construct the solution of JT]) we explore the relationship between the two 
evolution Eqs., fl) and ([7]). Specifically, we use the fact that any solution to ([TJ 
that satisfies ([2]) solves ([7]) as well, and that conversely, any solution to ([7]) that 
satisfies constraint (|2|) solves (fTJ). 

According to Lemma[2]the evolution of the constraint quantity under ([7|) is static. 
This implies that any solution to ([7]) corresponding to constraint-compatible data 
automatically satisfies for all t. Thus any solution to ([7]), (fT5|) obtained from 
constraint-compatible initial data solves ([I]), ([2]) as well. Existence is established. 
Similarly, multiple solutions to fTJ, ©; (H3 imply multiple solution solutions to ([7]), 
(|13p which is impossible in view of Theorem [3J Thus uniqueness is established. □ 

Remark. Theorem 0] states that the constrained evolution problem (jTJ) , |2]) is well- 
posed if the initial data is constraint-compatible and two boundary conditions are 
given in the form of (|13p (with freely specifiable gA )■ 

It is often advantageous from the numerical point of view to solve the free evolu- 
tion problem versus the fully constrained problem. We recall that only the evolution 
Eq. UJ is solved in the free evolution problem, while constraint ^ is monitored but 
not actively enforced on the solution. As is discussed in Section 2, to guarantee that 
(J2J) is satisfied for all times, the boundary data on variables of ([1]) must be given 
consistently with Eq. (j4|) that describes evolution of the constraint quantity. The 
following theorem gives an example of radiation-controlling boundary conditions 
that guarantee preservation of the constraint. 

Theorem 5. Let ctEl, |a| < 1, and the fields g and gA defined on the boundary 
be compatible at corners and satisfy 

(15) d t g = -d A g A , g(0) = (l + a)d t u n (0) + (l-a)d n u n (0), on dCl. 

Then there exists a unique solution to the free evolution problem {7]] satisfying the 
constraint compatible initial data (in the sense of Corollary Uty and the boundary 
conditions 

{d t u A + 2d[ n u A ]) = a(d t u A - 2d [n u A ]) + gA, 

(16) (d t u n + d„u n ) = -a(d t u n - d n u n ) + g. 

Moreover, the solution satisfies the constraint Eq. 0j. 



10 



ALEXANDER M. ALEKSEENKO 



Proof. To prove existence we recall that by Theorem|4]there exists a unique solution 
satisfying (p}, ([2|), and the first condition of (fT6l) . 

Next we expand C = d l Ui = d n u n + d A u A using normal and tangential deriva- 
tives and consider 

(17) d t C + d n C = d t d n u n + d t d A u A + d n d n u n + d n d A u A . 
Taking the normal component of ([I]) we have 

d\u n = d n d n u n + d A dAu n . 
Solving the last expression for d n d n u n and substituting into (|17p . we obtain 

d t C + d n C = d t d n u n + d t d A u A + d t d t u n - d A d A u n + d n d A u A . 
By commuting derivatives and rearranging terms in the last expression, we have 

(18) d t C + d n C = d t (d t u n + d n u n ) + d A (d t u A + 2d[ n u A] ). 
Similarly, 

(19) d t C - d n C = -d t (d t u n - d n u n ) + d A (d t u A - 2d [n u A] ). 

Let us show that the second condition of (fTB|) is a consequence of (JTJ) , (J2J), and 
the first condition. First of all since (J2J) is satisfied, both <9 t C + d n C = and 
dtC — d n C — 0. Multiplying (Tni|) by a and subtracting it from we obtain 

= d t [(d t u n + d n u n ) + a(d t u„ - d n u n j\ + d A [(d t u A + 2d[ n u A] )- a(d t u A - 2d [n u A] )}. 

In view of the first condition of (fTl)|) the latter can be rewritten as 

= d t [(d t u n + d n u n ) + a(d t u n - d n u n )\ + d A g A . 

Using (|15p we have (|16p by integration in time. Existence is proved. 

To prove uniqueness, it is sufficient to show that any solution to |T]), (Ti7i|) pre- 
serves the constraint, so that we can employ the second part of the proof of Theo- 
rem m 

By multiplying ^ with dtC, taking the integral over Q, and integrating by parts, 



(d;c)(d t c)dx = - / {d l c){d l d t c)dx+ / (d n c){d t c)da, 

i Jn Jon 

where is the outward pointing unit normal vector to the boundary. By commut- 
ing derivatives and reorganizing terms we rewrite the last equation as 

(20) d t [||d t c|| 2 L2(R3) + WdiCWl^} = \ Jj{d t c + d n cf - (d t c - d n cf]da. 

An argument similar to that of Corollary [1] can be used to show that if C satisfies 

(21) (d t C + d n C) = a(d t C - d n C) on Oil, 

with some constant \a\ < 1, the energy quantity in (|20p is not increasing. Thus the 
constraint quantity remains zero if it is zero initially. 

Indeed, (f2"Tj) follows from fig]) and (TT9")) by (fTg]) and dT5j). Uniqueness is proved. 

□ 

Remark. To establish Theorem [S] we use the assumption of the exact compatibility 
of data. Such assumption, however, is not practical in numerical applications where 
the equations are satisfied only approximately, and where one has to worry about 
propagation of small constraint violations. The next result generalizes Theorem [5] 
to the case of arbitrary data. 
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We rewrite Eq. |T]) as 

(22) d?u i =2d l d [l u i] +d i d l u l . 

We introduce the new variables 7Tj = dtUi, ipji — 2<9uiiji and p — d l ui, and decom- 
pose (J22) into 

d t Ui = 7Tj, 

d t TTi = d j ip oi + ditp, 

(23) dt<p = d l m. 

The resulting first order system is symmetric hyperbolic. Its characteristic speeds 
and fields in the direction rij are given by 

: if>AB, 

±1: — TT A ± IpnA, W ± = 7T„ ± if. 

Notice that the first two boundary conditions in (fl6|) remain the same for (l23l) and, 
in fact, take the familiar form of (jllj) with the substitution — aSij. The last 
condition, however, can not be implemented directly since d n u n is not among the 
variables tpij and p (but is closely related to ip). To formulate an equivalent con- 
dition by considering Eqs. (|23p at the face of the boundary we derive the following 
identity: 

(24) (d t <p ± 3 n p) = ±{d t -K n ± a„7T„) + d A {n A ± ip nA ). 
Let the differential operator L be defined by Eqs. (|23[) . i.e., 

L(iri,ipji,(p) = (d t TTi - d 3 ipji - df^j-i - 2dy-n; i \,dt<p - d l ni). 
As in Section 4, we define 

K L = {(TTj.Vji.fp) G L 2 (ft x [0,T]) | Lfruil/frtp) S L 2 (fi x [0,T]), 

|| (iri,1pji,ip) 1 1 « £ = 11(71-*, ^ji,¥?)|U 2 (Qx[0,T]) + ll i ( 7I " lI V'^,y)llL 2 (Ox[0,T])}- 

Theorem 6. Lei a € K, \a\ < 1, and the fields g and qa be defined on the boundary 
and compatible at corners and assume that g a meets the assumptions of Theorem^ 
Let, in addition, functions g and gA be such that (dtg + d A gA) G -ff 1 / 2 (9J7). Then 
there exists a unique solution (-Ki^ipji, ip) € Ti.^ to i23\) satisfying the boundary 
conditions 

w\ = aw A +.9.4, 

(25) (d t p + d n p) = a(d t p - d n ip) - (d t g + d A g A ). 
Moreover, the solution satisfies the energy estimate 

sup [2||7r(t)|| i2( n) + U(t)h 2 (n) + \\<p{t)\\ La{n) ] 

0<t<T 




+ [\\d l d t u l (0)\\L 2m + \\d l u l (0)h 2m ]) 
(26) +[2||7r(Q)|U a(n) + ||V(0)|| La(n) ], 

where the constants c, C\, c 2 are independent ofni, 'ipji, p. 
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Proof. By differentiating the last of Eqs. (|23|) with respect to time and substituting 
the second in the result, we obtain 

(27) %<p = d l d^. 

The initial data for §Ff\ is given by (p(0) = d l ui(0) and d t (p(0) = d l d t ui(Q). By 
employing the first order reduction, it is straightforward to establish existence and 
uniqueness of the solution. Indeed, by introducing d t ip — £ and Q — diip, we obtain 
the first order symmetric hyperbolic decomposition of (|27[). Notice that, when 
written in terms of £ and Q, f| 2 5 j) represents the maximally nonnegative boundary 
conditions (which are in general non- homogeneous). To establish the existence, it 
is sufficient to notice that the assumption {d t g + d A gA) € if 1 / 2 (90) guarantees 
the existence of functions h and hi such that £ = £ — h, Q = ~ hi satisfy 
the homogeneous conditions. Indeed, one could take h to be the H 1 extension 
of (d t g + d A gA)j c l 1 and set h\ — dip, where p solves d l dip = 0, {d/dn)p\eo, — 
{d t g + d A gA)/2- By completing the argument as in Section 4, we arrive at the 
estimate 

sup [\\dMt)\\ L2(a) + \\diip(t)\\ L2{n) ] < ci / \\dtg + d A g A \\ H i/2 {dn) dt 

0<t<T JO 

(28) +[11^(0)11^(^ + 11^(0)11^(0)]. 

Finally, by repeating the constructions of Theorems[2]and[lJ with the only difference 
being that di<p is now the forcing term in the second equation of (|23[) . we establish 
the existence result and the estimate (|26l) . The term ||v(*)||wn) ^ s added to (|26]) 
by considering the inequality 

^IMIi 2( o) = mw) < + M 2 

and using ([28]). □ 



Remark. It seems to be more natural to seek a solution to the problem (|23p in 
the space = {(iri,ipji,<p) € L 2 (fl) \ L(in,ipji,tp) G _ff~ 1 (fi)} (see [3S]). This 
is possible if one extends the results of 25J to weaker spaces and constructs the 
solution ip G Li2(£l), d l di(p £ H~ 2 (Q) to (|27p. A solution in this form does not 
require the "extra" smoothness assumption {dtg + d A gA) € H 1 ^ 2 {di}). Instead, it 
requires the milder assumption {d t g + d A gA) € H~ 1 / 2 {dVl) which is the consequence 
of the other assumptions of Theorem[6] It seems likely that one could use the general 
scheme of [25^ to prove the existence without the extra smoothness assumption. The 
proof, however, is rather technically involved. We, therefore, do not attempt it in 
this work. 

Theorem [6] is used to formulate a well-posedness result for the free evolution 
problem |Q}, (fT6| . 

Theorem 7. Let a G K, \a\ < 1, and the fields g and gA satisfy the assumptions 
of Theorem^ Let, in addition, the initial and boundary data be compatible, i.e., 



(29) g{0) = {l + a)d t u n {0) + {l-a)d n u n {0), on dQ. 
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Then a unique solution m G L2(0,), duUj] G ^(fi), d L ui G £2(^2) exists to i22\) 
satisfying boundary conditions 116\) . Moreover, the solution satisfies the estimate 

sup [2\\d t Ui\\ L2(n) + ||%w 2 ]||l 2 (q) + ll^u ; || L2(a) ] 




+ [\\d l dM0)\\ L2m + ||dV0)IU 2 (n)]) + [2||ft«i(0)|U a(n) + ||%^](0)|U 2( n)]- 

Proof. Existence of the solution follows from the first order reduction and Theo- 
rem [6] It should be noted, however, that the first order reduction ([23]) involves the 
constraint Cji — ipji — 2duU{\ = 0. Thus a solution to (f23|) only yields a solution to 
([22)1 if this constraint is satisfied. However, this poses no problems since (|23|) yields 
that dtCji = 0. Since Cji — for the initial data Ui(Q), if>ji(0), the latter implies 
that Cji = is automatically satisfied for future times. Thus the component Ui of 
the solution to (12T?)) gives a solution to (f2"2")) . Finally, one has to check that second 
part of (TTS|) is satisfied by Ui. In view of (f2T))) . the latter follows from identities (f2"4")l 
and conditions (|25p by integration in time. For the uniqueness proof, it is suffi- 
cient to notice that any solution to p2|) . (fl6|) solves (|23|) in the sense of the first 
order reduction described above. Moreover, the variable <p satisfies (|25j) in view of 
identities (|24[). Then the uniqueness follows from Theorem [6l □ 



6. The full gradient estimate for the static constraint evolution 

SYSTEM 

Energy estimates of the type (fl4)) are often considered to be a disadvantage 
for the numerical solution of (JTJ) , because they do not control all components of 
the gradient. The main objection to (fT4|) comes from the fact that if lower order 
(nonlinear) terms are present in the equations that include components of the gra- 
dient not controlled by ([HI) , coupling to such not-controlled terms may threaten 
the stability of the entire system. 

Next we will show how energy estimate (|14[) can be improved by employing 
the constraint preservation property. It has to be emphasized, however, that the 
existence and uniqueness results above suggest that energy estimate (|14j) is accurate 
for the problem and may not be improved unless one further restricts the class 
of solutions by specifying special boundary or initial data, or making additional 
assumptions about the equations (e.g., the gauge choice of [27] for the fat Maxwell's 
system). 

The following energy argument can be applied to a static constraint evolution 
system to provide additional energy estimates. Rewrite J7]) in the divergence form, 

dim = d J (djUi - Sjid p u p ). 

Contracting the last expression with dtu* and integrating by parts, we have 

7;dt\\d t u\\ 2 L (n) =- { {d-jUi- 5j l d p u p ){d : >dtu l )dx+ [ (d n Ui - n l d p u p )(d t u t )do-. 
1 ' hi Jon 
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Recalling that C = d p u p we expand the identity by adding and subtracting terms 
in the first integral: 

;9tllWi J (o) = - / {djui- SjiCXdidtv? - 5 ji d t C)dx 
Jn 

{djUi - S Jl C)(5 jz d t C)dx + / (<9„w 4 - nid p u p )(d t u l )da. 
n Jan 

The latter becomes after rearrangement 

^[||fle«lli a( n) + W d i u i " SjiC\\l a(p) ] = - J(djUi - 5 J% C)(5^d t C)dx 



(30) + / (d n ui- ni d p u p )(d t u l )da. 

Jan 

Let us for a moment assume that a solution to Eq. {7J) is defined in the entire M 3 
and is integrable. Then the boundary integral in the last expression drops out. 
Also, if the solution Ui is obtained from the constraint compatible initial data, then 
the first integral at the right side as well as all other terms in C drop out due to the 
trivial constraint preservation (see Lemma [5]). As a result, we obtain an estimate 
on all components of the solution's gradient. 

If the initial data is not constraint compatible, e.g., is perturbed by roundoff 
error, and one is dealing with a bounded domain Q, C K 3 one can rewrite (f30|) as 



dt[\\d t Ui\\l 2m + \\d jUi - 5 3l C\\i 2{n) - 2||C||£ 2(n) ] = 2 / (d nUi - ni d p u p )(d t u l )da. 

Jan 

Due to the static constraint evolution, C = d t C(0)t + C(0). In other words, ||C|| 
is completely determined by the initial data. Thus the left side gives the norm 
of the gradient of the solution up to some integral terms depending on the initial 
values Ui(0), d t Ui(0). We define e = ||<9 t %||£ 2(u) + \\djiH - <5jiC||! 2(n) - 2\\C\\l 2{n) 
noting that the right side of the identity contains a boundary integral and is thus 
dependent on the particular boundary conditions. 

Rewriting the integrand in terms of normal and tangential components and sim- 
plifying, we have 

(3i) d t e = 2 [ (d n u A d t u A + d A u A d tUn ). 

Jan 

It may look as if the boundary integral requires conditions on all three components 
of Ui but this is not so. In fact, the energy e is conserved if we require the solution 
to satisfy 



(32) u A 



= 0. 

an 



For, the integrand on the right side vanishes. Therefore the full gradient norm of 
the solution to ([7|), (|32|) grows at a rate not faster than t 2 . 

Another interesting example is given by the boundary conditions: 

(33) (d n u A - d A u n )\on = 0. 

Notice that taking the normal component of ([7]) at the boundary, we obtain 

d 2 u n = d n d n u n + d A d A u n . 
Due to the static constraint evolution, d l Ui = dtC(0)t + C(0), which implies that 
d n d n u n = -8 n d A u A + d n (d t C(0)t + C(0)). 
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Combining the last two equations and using (133p . 

d 2 t u n = d A d A u n - d n d A u A + d n (d t c(o)t + c(o)) - d n (d t c(o)t + c(o)). 

In other words, we have showed that the component u n at the boundary is deter- 
mined by the initial data. If the initial data is such that (a) u n \an = dtu n \ga = 0, 
and (b) the data is constraint compatible (so C = 0), then the last identity implies 
that u n \an = 0. Together with (J33J), this implies that the boundary integral in ([3~1|) 
vanishes. (Notice that assumption (a) overlaps with the well-posedness result of 
[27] for the fat Maxwell's system.) 

Remark. A similar full gradient estimate can be derived from the identity which is 
true for all H 1 vector functions (not necessarily solutions of ([7])): 

ll% u »]lli 2 (n) + ll^illwm + / (djUin 1 )^ = Wd^Wl (n) + / {d'u^ri'uj. 
' Jan Jan 

Using this identity, a full gradient estimate can be developed for the radiative 
condition (|13[) . The estimate does not need (a), but requires a special assumption, 
namely, the half plane reduction for the problem. We investigate this further in the 
next section. 

7. Comparison stability analysis by Laplace-Fourier transform 

In this section we make a comparison of the data derived by the static constraint 
evolution reduction to the constraint-preserving boundary conditions obtained by 
the standard methods using the subordinate constraint evolution equation (cf., 
[U [21 [U [TBI [El Ell [27] ) . In [27] a general class of constraint-preserving boundary 
conditions is derived for a symmetric hyperbolic reduction of the fat Maxwell's 
system and the well-posedness is proved for a subclass of these conditions (which, 
in fact, overlaps with (Til)]) in two equations). We will study (fT5|) and the conditions 
of [27] to see if any significant difference can be observed in the Lapalace-Fourier 
analysis. 

From ([20]) it is apparent that as soon as either of the three conditions is satisfied, 

(34) (7 = 0, d n C = 0, or d t C + d n C = 0, on dn, 

the constraint quantity is preserved. 

First, we consider the classical examples of Dirichlet-Neumann constraint-preserving 
boundary conditions (e.g., [TJ [2l [HI [32] ) . Let g A be a sufficiently smooth com- 
patible vector field tangential to the boundary. It can be easily verified that the 
following boundary conditions imply C\gn = 0: 

(35) u A = g Al d n u n = -d A g A . 

To construct the boundary condition that enforces d n C — 0, recall that 

d n C = d n (d n u n + d A u A ) = 0. 

By trading the second normal derivative of u n from ([T]) for the temporal and tan- 
gential derivatives, one obtains 

(36) d 2 t u n - d A d A u n + d n d A u A = on 90. 

One can either impose (|36| as a (second order) differential boundary condition, or 
to evolve it along the boundary using the initial data and compatibility conditions 
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at corners to produce the value of u n at the boundary. In this case, the constraint- 
preserving boundary condition takes the form 

(37) d n u A = h A , u n = u n , 

where h A is a given tangential vector field and u n is the solution of (136p . 

Finally, let us construct the examples of radiative conditions similar to one that 
can be found in [57] (in which the analog of (|3"T|) is also considered). Substituting 
C = d n u n + d A u A into the third condition of (|34[) , commuting derivatives, and 
trading the second normal derivative for the second temporal derivatives in p]). we 
have 

(38) d t (d t u„ + d n u n ) - d A d A u n + d A (d t u A + d n u A ) = 0. 

As in [27] we treat this equation as an evolution equation for dtu n + d n u n . It follows 
that the following conditions are constraint-preserving: 

(39) d t u A + d n u A = f A , d t u n + d n u n = /; on dfl. 

Here / = d t u n + d n u n is obtained by solving ([38]) . 

Unfortunately, the equation for the forcing term / does not decouple from the 
system. We will still include (f3"9"|) in the analysis since it can be implemented 
numerically. 

We now consider the half space x € R 3 , x% > and write the first order reduction 
of Eq. Q: 

d t u.i = 7r 4 ; ip ki = d k Ui; 

(40) d t ipji = djWi. 



Boundary conditions (|T6[) . (|35|) . (I37|) . (I39|) are imposed at x\ — 0. In terms of the 
first order variables these boundary conditions read (here g A , <?, h Al h, and f A are 
free functions) 

1 + a 

(41) 0ir x -tp xx =g, tta - 0{<PxA - <fAx) = 9A, #=:; , \a\ < 1, 



a 



(42) <p xx =h, n A = h Al 

(43) ir x = h, (p xA = h A , 

(44) d t {iT x - ifixx) - d A {f A + ipAx) = 0, Tr A -<p xA = f A - 

Notice that the right sides g A , g, h Al and h are no longer required to be compat- 
ible with the constraint. This reflects the situation common in numerical calcu- 
lations when the constraint-preserving boundary conditions are perturbed by the 
discretization errors. 

We seek a solution of (|4T)]) as a superposition of modes 

(45) ( n ) = ( ^\ )exp(st + iu A x A ), 

where s £ C, ^R(s) > 0, u> A is a real two-vector, and fields TTi{x\), (pji{x\) are 
differentiable. By inserting (|4"5"|) into the last two equations of ([40]) . we obtain the 
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system of ordinary differential equations: 

stt 1 {x 1 ) = d x if> xi (xi) + Ujj A y> Ai {x\); 
s(p xl (xi) = d x ix l {xi); 

(46) s(f>Ai{x\) = UVATTi{xi). 



Using the fact that 5R(s) > 0, we solve the last equation of for (pAi- Substituting 
the obtained expression into the first equation, we have 



(47) d x (p xl {x x ) = ([s 2 + \lo\ 2 ]/s)^{ Xi ). 



The system's matrix has eigenvalues fj/^ = ±\J s 2 + \u)\ 2 (the branch of the square 
root is chosen such that 3?(^ + ) > 0), and eigenvectors r~ 1 (svi,±^s 2 + \uj\ 2 Vi). 
Here Vi is a unit three- vector, r > is the normalization factor r 2 = |s 2 + |w| 2 | + |s| 2 , 
where \s\ = ^/~sE. 

The most general solution of (|4"7) that is bounded at infinity takes the form 



(48) ( )=J2^eM-V^+W^i), 



(f>xi{xi) 



3 



i/=l 



where w v = r 1 (s(v v ) i , -^/s 2 + |w| 2 (v„)i), and (u 1 ) i = (1,0,0), (u 2 )i = (0,1,0), 
( w 3)i = (0,0,1). Notice that in (|48f only the components corresponding to the 
eigenvalues with negative real part are taken. 

The coefficients a v are to be determined from the boundary conditions (|4T|) (|44l . 
By taking the Laplace transform with respect to time and the Fourier transform 
in variables x a = (^2 ; x^), and using the last equation of (|46[) to eliminate (p A %, we 
rewrite boundary conditions (|4Tj) - (|44"]) as 



6n x (0) - (p xx (0) = g, n A (0) - 0(tp xA (0) - — tt x (0)) - g A , 6 = \a\ < 1, 

s 1 — a 

fixxfi) = h, tta(0) = h A , 
Tr x (0) = h, (f xA {0) = h A , 

(49) s2 + ^ 2 n x (0)- S £ xx {0) = iLu A f A , n A (0) - ? xA (0) = h- 
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In each of the four cases the coefficients a v can be uniquely determined, correspond- 
ingly, by 
(50) 




Once the coefficients a v are found, the solution to ((40]) is obtained by substituting 
(TO into ((45]). 

The necessary condition for stability is that ((40]) must admit only the trivial 
solution satisfying the homogeneous boundary conditions [T31 H71 |T5] . Otherwise, 
with each solution (j4"5)) . the function 



= ~ > \ exuiast + lacuAX I 
V <PJi J V fjiiaxi) J 

is also a solution, for any a > 0. Since we can choose the parameter a to be arbitrar- 
ily large, the solution of (|40[) can not be bounded by any exponential function [T7] . 
The corresponding initial-boundary value problem is therefore strongly unstable. 

It can be seen from (|50p - (j5"5)) that if g A = h A = /a = 0, g = h = 0, then the 
coefficients a v = 0. Thus (jTi(xi), (p X i{x\)) = 0, and the necessary condition for 
stability is satisfied for each of the four sets of boundary conditions. 

In the case of vanishing initial data, well-poscdncss of the initial-boundary value 
problem for system ((4^]) follows if (]4"5]) is uniformly bounded by the boundary data. 
This is also known as the Kreiss condition [T31 [T71 [TB] . Specifically, if a constant 
K > independent of s and to a can be found such that (here || • || is the pointwise 
vector norm) 

(54) \\(M0),^(0))\\<K\\~g,9A\\ (< K\\h,h A \\, < K\\f A \\), 

then the unique solution to (|40p can be constructed by the inverse Laplace-Fourier 
transform, and can also be estimated by the boundary data [T3l [17] . 

It can be noticed that Kreiss condition is not satisfied by (|42[) and ((43]) . Indeed, 
by fixing uja and letting s — > 0, one can use ((BT]) . ((52]) to verify that the coefficients 
av can not be bounded uniformly in s. Also, the coefficients can not be uniformly 
bounded in u>a either. 

To discuss (j4"T]) and (j4"4"]) we will need the following lemma from [27] 

Lemma 5. Let P > 0, A, B £ R, (^4, B) ^ (0, 0), and consider the function 

V : MO > 0} -> c, ii>(0 = M-BVe + P 2 , 

where we choose the branch of the square root such that 5?(a/C 2 + P 2 ) > for 

»(c) > o. 
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Then, ip has zeros if and only if A > B > or A < B < 0. Furthermore, \ip\ is 
uniformly bounded away from zero if and only if AB < 0. 

From Lemma [5] it is apparent that as long as \a\ < 1, both s 2 + \ia\ 2 + 6s 
and Q\J s 2 + \lo\ 2 + s are uniformly bounded away from zero. This implies that 
coefficients a v in (|50| remain bounded if s — > 0, 5R(s) > 0. At the same time, when 
|s| -> oo, 

r _ y^TH 2 ! + \s\ 2 _ + H 2 /s 2 | + 1 2 

Similarly, when |c<j| 2 — > oo, the above expression converges to 1. Therefore, we 
conclude that <ti is uniformly bounded by g. By a similar argument, 02,3 are also 
uniformly bounded by (72,3 and g since the coefficient in braces can be treated 
similar to that of 01, and also, |w|(-\/s 2 + |cj| 2 + Os)^ 1 — > 1 as |w| — > 00. Thus the 
Kreiss condition is satisfied for (|4T|) . 

The condition (pH)) does not meet the Kreiss condition since o\ in (|55)) depends on 
1/ ^Js 2 + \uj\ 2 , which can be made large by choosing s = |e| + £ — > 0. However, 
the coefficients remain bounded in the high frequency regime, when \lu\ — ► 00. 

Remark. The necessary condition for stability is equivalent to the non-vanishing 
determinant of the corresponding matrix for the coefficients ov (cf., [13l[T7l[T8ll28| ). 
At the same time, for the Kreiss condition the explicit formulae for the coefficients 
a v are more informative. In particular, in ()50[1 . the matrix appears to be lower- 
triangular and thus its determinant does not include off-diagonal terms. As the 
result, terms {9iu>2^){ \J s 2 + \lo\ 2 + Os)^ 1 in the expressions for cr 2 ,3 are overlooked 
in the determinant condition. However, they clearly depend on loa- 



8. Weak instabilities 

A framework to quantify instabilities that grow polynomially in time was for- 
mulated in [28j . Such instabilities may be triggered by small perturbations of the 
initial data, and, for certain types of boundary conditions, can grow faster in time 
than a polynomial of any degree. We complement the analysis of Section 7 with 
this technique to see if non-vanishing initial data can produce a milder instability, 
similar to the one described in [2"T] . 

Following [28], we perform the Fourier transform of the solution in tangential 
variables xa and look for a solution in the form 



(55) 



* =exp(zu W A )]T : {u) ) 1 



where u S K, to > 0, loa is a real two- vector, (fr| , (p^ ') 7^ 0, and p > 1. 
Substituting {55]) into (jlu]) , and grouping terms by powers of i, we have 

t p : 



(56) = i^'twn) + w A ^}{uxx), 
= d x ni p) (tu Xl ), 

(57) = *<W p) (^i); 
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r,v = 0,...,p-l: 

(58) = d x ^(cj Xl ) + w^M, 

(59) ^ 1) (u;x l )=d x n^(cjx 1 ), 

(60) = iuj a tx\ v \ujxi). 

Obviously, 7rj = as long as ^ 0. We let ^ = p~ 1, differentiate ([59]) in n, 



and substitute the result in (ITO)) . Similarly, we substitute ([6"0")) in ([56")) to obtain 

= %^ 1 \ux x )-\w A \H'r l \u>x x ). 
The most general solution that is bounded at infinity is 

3 

^'''(wii) = ^2 o-i P_1) w Q exp(-|w j4 |a;i), 

a=l 

where ui = (1, 0, 0), v 2 = (0, 1, 0), and v 3 = (0, 0, 1). 

The coefficients <ri P 1 ^ are to be determined from the boundary conditions in 



view of the following relations that are consequences of (|59j) and (|60|) 

(61) tp%(uxx) = -|wa|^ (p_1) (wxi), ^(wh) = iu; A ^i (i^xi). 

One can check that <p^(ujxi) and (o>a;i) obtained from (f6Tj) satisfy (|56[) auto- 
matically. 

Similarly, for i/ = 2, . . . ,p, one finds 

~(v)( \ q2~(i / — 2)/ \ I |2-(^— 2)/ \ 

which can be recursively solved. As in the previous case, the general formulas 
for tp^) and <p%) can be obtained from (f59"|) and (frJO")) . and the coefficients are 
determined from the boundary conditions. 

Consider the case of boundary conditions (|4Tj) . By taking the Fourier transform 
in xa of gA and g, and decomposing the result in the Taylor series, we have 

(62) <Wr(">(0) - £g>(0) = §M, t#>(0) - 0(^j(O) - ^(0)) = §M 
where f = 0, 1, . . . ,p, and 9 = j^, \a\ < 1. 

Let functions (7,4 and g be such that = 0, <?( p ' = 0. If a nontrivial solution 
([55)1 (with some (p^) 7^ 0) can be found satisfying boundary conditions ([52]). then 
the initial-boundary value problem is ill-posed. Indeed, in this case, the solution 
at times t > is 0(uj p ) whereas the initial data is 0(u>°) and the boundary data is 
O^p^ 1 ). Therefore, the solution can not be estimated in terms of the initial data 
as u> — > 00, and the problem is ill-posed 28J. 

Using 7t£ = and (fBTj) we have 

f63l ct (p_1) - — 5W a [p - 1] - 1 5 (p) - if^iigW 

\ 00 / "1 — 1 iff ' "2,3 — /)i ii/2,3 1 12" 

This guarantees that 7^. and thus ^(wii) and ip^(ioxi) are zero as 

long as = and = 0. Therefore, (|4Tj) is polynomially stable. 
At the same time boundary conditions (|42|) yield, 

^(Q) = ftW, ffi l ' ) (0) = ^ ) I i/ = 0,l,...,p. 
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It can be seen that if = and = 0, but 7^ 0, the second condi- 

tion implies that ^(luxi) = ^ exp(—\ujA\xi) i= 0, and, in view of (|61[) . 
(pf\(ujxi) = — IwaI^-^ -1 ^ exp(— IwaIxi) ^ 0. Thus condition (|42j) admits ill-posed 
polynomially growing modes (unless Ha = 0, h = 0). 

Similarly, (|43|) is ill-posed, since it implies that if Z^ 23-1 ) 7^ 0, then (luxi) = 
-IvAlhb-V exp(-\w A \xi) ^0. 

Finally, consider boundary conditions (I44| . Using the Fourier transformation in 
variables x a and substituting decomposition (|55|) . the boundary condition yields 

(64) *%\0)-<P%(P)=ft\ i/ = 0,l,...,p, 

(65) ^(/W + ^(0))=0, 

(66) ttM(O) - 0M(O) - ^(/j^ + ^" 1} (0)) = 0, «/ = 1,2, . . . ,p. 

Since = 0, from §Q it follows that (p^ ] A = if /\ p) = 0. Therefore, from ([ST]). 
^A = and = 0. It remains to check that (pf} — 0, and (f>A x = 0. First 
of all u> A tf { Al = due to (65]). This, in view of ([61]), implies that ni P ~ 1] = and 
(f>xx — 0. By ([6"T]) this yields </5^2 = as well. Thus boundary condition (|44|) is 
stable in the sense that no polynomially growing modes that can not be estimated 
by the data is generated. 
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